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Abstract. We have obtained optical {BVR, 
[Oiii]A5007Aand Ha), near infrared [JHK) imaging and 
long-slit optical spectroscopy for the Blue Compact Dwarf 
galaxy Mrk 86 (NGC 2537). In this paper, the first of two, 
we present optical-near-infrared colors and emission-line 
fluxes for the currently star-forming regions, intemediate 
aged starburst and underlying stellar population. We 
also describe the evolutionary synthesis models used in 
Paper II. 

The R and Ha luminosity distributions of the galaxy 
star-forming regions show maxima at Mj^^— 9.5™ and 
-^'Ha=10^'''^ ergs~^. The underlying stellar population 
shows an exponential surface brigthness profile with cen- 
tral value, /iE, 0=21.5 mag arcsec"^, and scale, a=0.88kpc, 
both measured in the i?-band image. In the galaxy outer 
regions, dominated by this component, no significant color 
gradients are observed. 

Finally, a complete set of evolutionary synthesis mod- 
els have been developed, covering a wide range in metallic- 
ity, 1/50Zq<Z<2Zq, and burst strength, I-IO"". These 
models include nebular continuum and recombination and 
forbidden-line emission. 



that cither these objects are truly young systems or they 
have an intermittent star formation history with short 
intense star-forming episodes followed by long quiescent 
phases. Although this ques tion remained unanswere d dur- 
ing decades (Th uan |1983| ; Campbell & Terlevich |1984 ; 
Loose & Thuan |1985 ), nowadays, most of the studies 



on dwarf star-forming systems, including Blue Compact 
Dwarf galaxies (BCD hereafter), have revealed the exis- 
tence of an evolved underlying population (Kunth et al. 
1988t Hoff man et al. |1990| ; Papaderos et al. |l996b[ Dou- 



blier 1998; Norton & Salzer 1997) 



Therefore, the understanding of the mechanism, or 
mechanisms, governing this star formation regulation is 
the stepping stone of the dwarf star- forming galaxies evo- 
lution. Moreover, this mechanism constitutes the missing 
link between the diff erent dwarf star-forming gal axies (Silk 
et al. 119871 Burkert |1989[ Drinkwater & Hardy |99ll Pa- 



paderos et al. 1996b) 



Key viTords: individual: Mrk ; 
clusters, stellar content 



There are two different approaches to address these 
questions. On the one hand, statistical analysis of a sample 
of these objects will allow the study of the relationships be- 
tween fundamental parameters and properties of the star- 
for ming dwarf ga laxies (Marlowe et al. 1995 ; Papaderos et 
galaxies: irregular, star al. 1996a , 19961). On the other hand, the detailed analy- 



1. Introduction 

The high star formation rate (SFR hereafter) and low neu- 
tral gas content deduced for dwarf star-forming galaxies 
imply consumption time-scales of about 10^ yr (Fanelli 
et al. 119881; Thuan & Martin 119811), much shorter than 



the age of the Universe. Searle et al. (1972) suggested 
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sis of individual low-redshift objects with multiple regions 
of star formation is fundamental to reconstruct their star 
formation histories. In particular, we can obtain a bet- 
ter understanding of possible effects of merging or some 
internal process as self-propagation that could originate 
the spreading of star formation, and the effects of these 
star-forming events on the future star formation. 

Among dwarf galaxies. Blue Compact Dwarf galaxies 
conform the subpopulation where the star-forming events 
are most violent. Blue Compact Dwarfs are low-luminosity 
(Ms > — 18™ )|^ galaxies with compact sizes whose spectra 
are similar to those of low metallicity H ii regions (Searle 
fc Sa rgent |1972| ; Kunth & Sargent |1986| ; Thuan & Martin 
1981). Their spectra are characterized by emission lines 
over a blue continuum which implies the existence of a 



For _ffo=75kms-i Mpc"^ (Thuan & Martin 1981) 
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large fraction of OB stars and an intense star-forming ac- 
tivity. 

The Blue Com pact Galax y Mrk 86=NGC 2537 (Shap- 
ley fc Ame s 1932| ; Markarian 1969), also known as Arp 6 
(Arp 1966 ), constitutes an excellent laboratory to test the 
BCD star formation history since its star-forming regions 
populate all the galaxy. This object is the propotype of the 
iE galaxies, the most important BCD galaxies subclass, 
that co nform 70 per cent of the BCD galaxies (Thuan & 
Martin |l98lt Thuan |l99l| ). 

Optical imaging (both broad- and narrow-band) and 
near infrared images have been obtained to map the differ- 
ent regions and get insight into the star formation history 
of every region. In addition, we have obtained very high 
quality spectra with spatial resolution. The spectra of 22 
star-forming knots and of free interim were taken with 10 
distinct slit positions and orientations and moderate and 
high spectral resolution. 

After briefly introduce the previous results on Mrk 86 
in Sect. ^, we will describe the observations and data 
reduction in Sect. |^. The data analysis is presented in 
We describe the evolutionary synthesis models in 
In Sect. 1^ we give some results. Finally, our con- 
clusions are presented in Sect. ^ 

In Gil de Paz et al. (2000a[ Paper II hereafter) we 



Table 1. Collection of the published parameters of Mrk 86 




derive the physical properties of the star-forming regions 
and stellar populations of Mrk 86 comparing our photo- 
metric data with the evolutionary synthesis models pre- 
sented in this paper. Electron densities and temperatures 
and chemical abundances of the ionized gas in the galaxy 
star-forming regions will be also obtained. Finally, a global 
interpretation of the galaxy will be given. 



2. Mrk 86 

Mrk 86 (a(2000)=08'n3"14.56^ (5(2000)=-l-45°59'30.2")Q 
is a well-known object among BCD galaxies (Loose & 
Thuan |1985|; Sage et al. |1992| see Table |l]). In the Loose & 
Thuan (|l98§) BCD classification, Mrk 86 belongs to the iE 
class: smooth elliptical Low Surface Brightness (LSB) un- 
derlying stellar component on which several knots of star 
formation are superimposed. The present star formation 
activity is spread out in several clumps over a relatively 
large fraction of its entire surface. In fact, Mrk 86 shows 
the highest ratio of the area covered by star-forming re- 
gions to the total projected area o f the galaxy within the 
sample of Papaderos et al. ( 1996a ). 

To this moment Mrk 86 has not been the subject of 
a single and proof study, inasmuch the large number of 
published data on it. A collection of the published data 
IS given m Table |. A more detailed description of the 
previous observations is given by Gil de Paz et al. ( pJ99| ; 
GZG hereafter). 



Parameter 


Data 


Reference 


V© 


447 km s^' 


(1) 


vlg 


522kms-i 


(1) 


VLG 


460 km s"^ 


(2) 


Distance 


6.9 Mpc 


(3) 


B 


12.8 


(4) 


f(1510A) 


0.87X 10"" ergm"^ s"^ A"^ 


(5) 


S(2.8cm) 


7±2 mjy 


(6) 


S(6.3cm) 


18±4mJy 


(7) 


S(1.2cm) 


ll±4mJy 


(7) 


f(12/im) 


0.25 Jy 


(8) 


f(25/im) 


0.42 Jy 


(8) 


f(60/im) 


3.15 Jy 


(8) 


f(lOOAim) 


6.26 Jy 


(8) 




0.35 10** Lq 


(8) 


J 


13.25 


(9) 


H 


12.57 


(9) 


K 


12.37 


(9) 


L 


12.63 


(9) 


M(Hi) 


2.4 10** M0 


(10) 


Mt 


5.810** Mq 


(10) 


f(Hl 21cm) 


2.1Jy kms"^ 


(11) 


FWHMhi 


93kms"^ 


(11) 


Ico 


0.83±0.12Kkms-^ 


(12) 



(1) Bottine lli et al. ( |l990| ); (2) Rood & Dick cl (119761) ; (3) Sha- 
rina et al. Jl999| ); (4) Dultzin- Hacya n et al. ( |1990D ; (5) Fanel li 

_( [l99l|) ; (7) Kle in et al. ( |l984| ); 
Thuan (|l983| , measured 
Verter 



et al. ( [1988| ); (6) Klein et al 
(8) Thronson & Telesco (|l98^)1 
with 7 " 8 apertures); (10 ) Thu an 
(|l985|); (12) Sage et al. (199^). 



(9) 

& Martin (|l98l|); (11) 



^ Cente r of the outer isophotes as measured in the 7?-band 
(see Sect. 4.1) 



2.1. Distance 

Mrk 86 is a nearby star-forming galaxy. The heliocentric 
recession velocity is W0=447kms^^ (see recession veloci- 
ties referred to the Local Group in Table p. 

Using the mean blue magnitude of the three brightest 
blue stars Sharina et al. (1999) derived a distance for this 
object of 6.9 Mpc with an uncertainty of 20 per cent (M. 
E. Sharina, private communication). We can also estimate 
the distance to Mrk 86 if we assume that it is gravita- 
tionally bounded to the edge-on spiral galaxy UGC 4278 
(see Paper II). UGC 4278 has a corrected /-band abso- 
lute magnitude of -18.21±0.11 (see Giovanelh et al. [1997| ) 
with a logarithmic 21cm line width of 2.220±0.019 (in 
kms~^). Therefore, with respect to the Giovanelli et al. 
( 1997D TuUy-Fisher template, this galaxy shows an offset 
of 0.66™, which corresponds to a peculiar radial velocity 
of -253 km s-^ Thus, its recesion velocity in the CMB 
reference frame should be 968 km s^^ (wQ=565kms^^), 
which implies a minimum distance for Mrk 86 of 13 Mpc 
(for i/o=75kms~^ Mpc~^ and assuming that the differ- 
ence in heliocentric velocity is due to orbital motion) . The 
large difference between this value and that given by Sha- 



ina et al. ( 1999 ) suggests that these objects are probably 



A. Gil de Paz et al.: Mapping the star formation history of Mrk 86* 



3 



not physically bounded. Thus, since no other distance in- 
dicator is available for Mrk 86, we have used the distance 



Table 2. Journal of observations. 



measured by Sharina et al. (1999) throughout this work. 
This value leads a projected spatial scale of 33 pc arcsec^ ^ . 

3. Observations and reduction 

3.1. Optical imaging 

An optical _B-band image was obtained at the 1-m Ja- 
cobus Kapteyn Telescope (JKT) of Roque de los Mucha- 
chos observatory (La Palma, Spain) in 1997 November 
with a 24 ^m 1024 x 1024 pixels Tek#4 CCD (see Table |). 
An additional Johnson- image was taken at the 1.52-m 
Spanish telescope at EOCA (Calar Alto, Almerfa, Spain) 
in December 1993 with a Tek 1024x1024 CCD with pixel 
size of 19 /xm. This observation was broken up into five ex- 
posures with a total integration time of 2400 s. Finally, an 
i?-band image was obtained during service time in Novem- 
ber 1998 with the ING wide-field camera (ING-WFC) 
equipped with four EEV42 2kx4k pixels CCD detectors 
on the 2.5-m Isaac Newton Telescope (La Palma, Spain). 

Narrow-band images in the light of [Oiii]A5007A (Aq 
= 5012A, FWHM = 50A) and Ha (Aq = 6568A, FWHM 
= 95A) were obtained. In order to subtract the contin- 
uum, V and i?-band images were used, respectively. The 
[Oiii]A5007A image was secured for us on December 1993 
during service time with a 1280x1180 pixels EEV5 CCD 
attached to the 2.5-m Isaac Newton Telescope. The Ha 
image was taken during the same service time observa- 
tions that the i?-band image in November 1998 using the 
ING-WFC camera at the 2.5-m Isaac Newton Telescope 
(La Palma, Spain). 

High S/N ratio dome flats and exposures of the sky 
taken in twilight were obtained in any case. The standard 
procedure of bias removal, dark-current subtraction, and 
flat fielding using dome and sky fiat-field images was per- 
formed using the ESO image processing system MIDAS 
for the V image and IRAfQ for the B and i?-band im- 
ages. 

Atmospheric conditions were photometric during the 
observing runs. B and V images were flux calibrated ob- 
serving repeatedly during the nights a set of standard 
stars, taken from the lists given by Kent (1985) and Lan- 
dolt (1973). Finally, the i?-band image was flux calibrated 
using the radial surface brightness profile published by 
Papaderos et al. ( 1996a ). 

The [Oiii]A5007T and Ha images flux calibration was 
performed as follows (see also GZG). The continuum emis- 
sion was subtracted, using the V and R images, respec- 
tively. Then, we trimmed those regions covered by the slits 
of the 6(blue) and r(red) spectra (see Table 0). The b and 



Spectroscopic Observations' 





Exp. 


Slit 


Range 


Disp. 


Telescope 


time(s) 




(nm) 


(A/pix) 


CAHA2.2m 


3600 


#l,2,4,6fe 


330-580 


2.6 


CAHA2.2m 


1800 


#56 


330-580 


2.6 


CAHA2.2m 


3600 


#l,2,4,5,6r 


435-704 


2.6 


CAHA2.2m 


3600 


#3 


390-650 


2.6 


INT2.5m 


1800 


#7,8R 


637-677 


0.39 


INT2.5m 


900 


#9R 


637-677 


0.39 






Image Observations 






Exp. 


Filter 


Scale 


PSF 


Telescope 


time(s) 




("/pixel) 


(") 


JKTl.Om 


600 


B 


0'.'33 


I'.'O 


CAHA1.5m 


2400 


V 


0'.'33 


I'.'e 


INT2.5m 


900 


R 


0'.'333 


1'.'2 


INT2.5m 


900 


[Oiii]A5007 


0'.'57 


2'.' 5 


INT2.5m 


7200 


Ha 


0'.'333 


1'.'2 


KPN02.3m 


900 


J 


0'.'66 


I'.'S 


KPN02.3m 


360 


H 


0'.'66 


I'.'e 


KPN02.3m 


540 


K 


0'.'66 


1'.'7 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



See Fig. 1 of GZG for slit orientations and positions. 



r spectra were convolved with the transmission curves of 
the corresponding narrow-band filters. Then, we compared 
image counts and fluxes measured in the convolved spec- 
tra for several adjacent regions in order to determine the 
precise location of the slits. Finally, once the regions cov- 
ered by the slits were precisely defined, we corrected these 
calibration relations for the sensitivity of the filter at the 
corresponding wavelength and, in the case of the Ha flux, 
for the contribution of the [Nii]A6548A and [Nii]A6583A 
lines. The discrepancies obtained among different spectra 
were in all cases lower than 10 per cent. 

3.2. Near-infrared imaging 

Near-infrared images (nIR hereafter) of Mrk 86 in J (Aq — 
1.25/zm, FWHM = 0.30/im), H (Aq = l.Q^^iuv, FWHM = 
0.28Aim) and Ks (Aq = 2.15Aim, FWHA4 = Q.ii^im) were 
obtained on January 1998 with the Steward Observatory 
near-infrared camera equipped with a 256x256 NICMOSS 
detector attached to the 2.3-m Bok Telescope at Kitt Peak 
National Observatory (Arizona, USA). The observational 
proced ure closely follows th at of Ar agon-Salamanca et al. 
( |1993| ) and Gil de Paz et al. ( |2000bD . The total integration 
time of each subimage was broken up into background- 
limited sub-exposures to avoid saturation of the detector. 
Images of adjacent blank areas of the sky were alternated 
with the target frames for accurate flat-fielding. Compa- 
rable amounts of time were spent imaging the source and 
the sky to ensure adequate monitoring of the sky changes. 
Individual subexposures of the object were offset several 
arcseconds to improve the final result. 

The reduction process include: (1) bias and dark sub- 
traction of object and sky frames; (2) flat-fielding us- 
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Fig. 1. BVRJHK, Ha and [Oiii]A5007A images of Mrk 86. 



ing normalized sky frames; (3) sky subtraction using sky 
frames taken before and after each exposure; (4) bad pixel 
removal; (5) registering of the subimages using fractional 
pixel shifts and (6) median combining of all individual 
frames. The reduction was carried out using own IRAF 
procedures. 

The nIR images were calibrated o bservi ng standard 
stars from the list from Elias et al. ( 1982 ) during the 
nights at the same airmasses than the object. We have 



assumed a color independent correction between the Ks 
and K' bands {K' = KM ; Ap = 2.12/im, FWHM = 
0.34/xm; Wainscoat & Cowie 1992). In order to check the 



validity of this assumption, we convolved the Ks and K' 
filters response functions with Planck spectral distribu- 
tions at different temperatures in the range 3000-20000 K, 
obtaining the fluxes Fks and Fk'- The larger difference 
in 2.5x\og{FKs/ Fk') within this temperature range was 
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Fig. 2. Knot contours derived using the COBRA program superimposed on the Ha image. The thick-lined contour 
corresponds to the starburst component as described by Papaderos et al. (1996a, see Sect. |6.3| ). The crossed-circle 
position marks the center of the outer i?-band isophotes. 



0.025™, small enough to assume this correction to be in- 
dependent of the spectral energy distribution. 

Thus, once the fluxes were transformed to K'-hand 
fluxes, we converted them to the standard i^-band (Aq 
= 2.19^m, FWHM = 0.41^m) applying the s tandard 
correction given by Wainscoat & Cowie (1992), K' — 
K=0.22x{H - K). 

In Fig. |l] we show the Mrk 86 neighborhood in 
the optical-near-infrared bands studied, including 
the Ha and [Oiii]A5007A images. The astromet- 
ric calibration was performed using the program 
PlateAstrom (Garci'a-Dabo fc Gallego 199S|; see als o 



http : //www ■ ucm . es/info/Astrof /opera/opera . html ). 

The bright field star placed at the relative position 



(60"E,30"N), saturated in the V and i?-band images, was 
artificially removed. 



3.3. Spectroscopy 

Up to 14 optical long-slit spectra at 10 different slit posi- 
tions were obtained. The 11 medium-low resolution spec- 
tra (6 and r spectra; see Table |) were obtained with the 
Boiler & Chivens spectrograph at the Cassegrain focus 
of the 2.2-m telescope at the German-Spanish Calar Alto 
observatory (Almerfa, Spain) in January 1993. Althought 
special care was taken to place the slits into position, us- 
ing offsets from the bright field star at north of the galaxy, 
there were some lack of precision. The actual positions of 
the slits (see Fig. 1 of GZG) were determined a posteriori 
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Fig. 3. Surface brightness profiles in the BVRJHK 
bands. The fit to the exponential component of the R- 
band profile at distances larger than 70" is also drawn. 
Error bars represent ±lcr. The K-hand profile has been 
offset —0.5™ in order to avoid confusion with the _ff-band 
profile. 



comparing broad-band image spatial cuts with the spa- 
tial profiles of the spectra previously convolved with the 
corresponding filter transmission curves. The detector em- 
ployed was a 1024x 1024 Tek#6 CCD with a pixel size of 
24 /im. The 600grmm^^ grating chosen provided a spec- 
tral resolution of 6A in the light of Ha and a reciprocal 
dispersion of 2.6A/pixel, with a slit width of 2'.'65. The fi- 
nal spatial scale was l'.'43/pixel. The spectral coverage was 
around 2500A and the grating angle was selected to cover 
the blue region (~3300-5800A) and red domain (-4350- 
7045A) in two different exposures which overlapped on the 
H/3 region. The seeing was variable during the observing 
run with FWHM between l'.'5-2'.'5. The low airmasses at 
which these spectra were obtained (<1.2) guarantee that 
no significant loss of blue light due to atmospheric refrac- 
tion has ocurred. In addition, three high resolution spec- 
tra (7R,8R,9R) were obtained with the IDS spectrograph 
at the Isaac Newton Telescope (INT) of the Roque de 
los Muchachos Observatory (La Palma, Spain) in January 
1998 with a 1024x1024 Tek#3 24 /xm CCD (see Table |). 
The 1200R grating (1200grmm~^) was used with an slit 
width of 1", providing a spectral resolution of 0.9A in the 
light of Ha and a reciprocal dispersion of 0.39A/pixel. The 
spatial scale was of 0.33"/pixel. 



These spectra were reduced using the FIGARO (January 
1993) and IRAF (January 1998) software packages. After 
bias removal and flat-fielding, the frames were cleaned of 
cosmic rays. The sky was removed from each frame by 
subtracting a polynomial fitted to those regions free for 
object emission. Wavelength calibration was performed by 
using He-Ar lamps observed inmediately before and after 
the galaxy integration. The standard stars Hiltner 102 and 
Hiltner 600 were observed at different airmasses in order to 
correct for atmospheric extinction and to ensure absolute 
flux calibration. 

Finally, we requested an UV spectrum of Mrk 86 from 
the International Ultraviolet Explorer (lUE) Final Archive 
(see Fig. It was originally take n by A Uoin and Duflot 
in January 1983 (see Bonatto et al. 1999| ). The total expo- 
sure time of this spectrum, SWP18927, was 24000 s. It was 
obtained in low dispersion mode with the SWP camera. In 
the observing spectral range, 1150-1975A, the resolution 



power varies between 270-300 (Cassatella et al. 1985) 



4. Analysis 

4.I. Surface brightness profiles 

We have analyzed the surface brightness profile of Mrk 86 
in different bands {BVRJHK; see Fig. ^). These surface 
brightness profiles have been obtained using the IRAF 
task ELLIPSE. Background galaxies and foreground stars 
were interactively masked. We obtained the surface bright- 
ness profiles over the original images, excepting in the 
outer regions of the _R-band profile where a median fil- 
ter of 5x5 pixels (O'.'333/pixel) was appUed. Due to the 
contribution of about 60 high surface brightness regions, 
we only fitted the isophotcs with major axis larger than 
25". Using the mean of the isophotal centers measured be- 
tween 80 and 120" in the i?-band image we have estimated 
the galaxy coordinates given in Sect. I (see also Fig. |). 

4-2. Knot positions and sizes 

In order to derive the positions and sizes of the regions in 
the neighborhood of Mrk 86 we have developed a program 
called COBRA (see Appendix) . Briefly, this program selects 
the image section where the region of interest is placed. 
Then, the light profiles in both image axes are fitted us- 
ing a straight line, in order to account for the underlying 
emission, and two gaussians reproducing the knot emis- 
sion profile. The center of the knot is estimated as the 
maximum of this latter component. The positions derived 
for all these regions are given in Table ^ (columns 3 and 
4). They have been numbered in reverse declination or- 
der. Most of these regions were identified on the Johnson 
i?-band image, the deepest of those shown in Fig. || (For 
i?=25.6 the signal-to-noise ratio is 30). However, some re- 
gions showed intense Ha emission, but practically no op- 
tical or near-infrared continuum emission (see Table 



Table 3. Positions and sizes for the Mrk 
regions. 
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Table 3. (cont.) Positions and sizes for the Mrk 86 neigh- 
borhood regions. 



# 


#' 


RA(2000) 


DEC(2000) 




e 


CI. 


(1) 


(2) 




(3) 


(4) 


(5) (6) 


(7) 


(8) 


01 




08 


13:16.08 


+46:00:17.1 


0.71 PLR 


2 


B 


02 


— 


08 


13:16.33 


+46:00:12.1 


0.85 0.45 


1 


B 


03 


— 


08 


13:15.91 


+46:00:11.8 


0.74 0.17 


2 


B 


04 


01 


08 


13:14.72 


+46:00:04.1 


1.13 0.87 


2 


F 


05 


— 


08 


13:13.38 


+46:00:01.7 


0.91 0.56 


2 


B 


06 


02 


08 


13:12.36 


+45:59:57.6 


1.07 0.79 


2 


E 


07 


- 


08 


13:15.70 


+45:59:57.4 


0.83 0.41 


2 


E 


08 


03 


08 


13:16.18 


+45:59:56.0 


0.96 0.63 


2 


E 


ogt 


_ 


08 


13:15.28 


+45:59:54.1 


0.76 0.24 


2 


E 


lot 


_ 


08 


13:13.56 


+45:59:53.4 


0.68 PLR 


1 


E 


11 


04 


08 


13:12.20 


+45:59:53.0 


0.75 0.21 


2 


B 


12t 


_ 


08 


13:16.74 


+45:59:52.9 


1.06 0.78 


2 


E 


13 


05 


08 


13:14.36 


+45:59:52.0 


1.31 1.09 


1 


E 


14 


06 


08 


13:11.72 


+45:59:50.7 


1.36 1.15 


2 


S 


15 


08 


08 


13:13.85 


+45:59:50.7 


0.96 0.63 


2 


E 


16 


07 


08 


13:14.77 


+45:59:50.6 


1.72 1.56 


1 


S 


17 


— 


08 


13:09.23 


+45:59:49.5 


1.03 0.74 


2 


E 


18 


09 


08 


13:13.08 


+45:59:48.2 
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(1) Knot number (reverse DEC sorted). 

(2) Old knot number as given in GZG. 

(3) RA(J2000). 

(4) DEC(J2000). 

(5) Radius (in arcsec) at 1 e-folding. 

(6) Radius (in arcsec) corrected for the atmospheric seeing. 

(7) Folding (e'^, e^ or e^) where colors were measured. 



(8) Classification (see Sect. 4.3). 
PLR= Point-Like Region. 

+ Knot sizes measured on the Ha image (see Sect. 



Since more than thirty new objects are identified in this 
i?-band image relative to the Gunn-r image used in GZG, 
we have introduced a new notation for the knots. In this 
sense, the regions numbered as #9, #15, #16, #31 and 
#49 in GZG, now are, respectively, #18, #26, #27, #42 
and #70 (see Fig. | and Table |). 

The determination of the boundaries of these knots 
is not an obvious task. Some authors simply locate the 
position of the knots by using a weighted mean of the 
pixels surrounding a maximum of intensity and then they 
use large enough circular apertures to perform the pho- 
tometry. This simple approach assume circular geome- 
try for the knots, a prerequisite which is not always ful- 
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filled. An autonomous feature, commonly an H ii region, 
can be also defined as the region which is inside its out- 
ermost closed contour (Petrosian et al. 



LO 



1997) or select- " 



ing all pixel that are connected above a given threshold 
using predetermined background and noise properties of 
the frame (Mazzarella et al. 1993| ). All these methods, 



rewieved in Fuentes-Masip (1997), are difficult to apply in 
overcrowded fields or when there is a very intense and vari- 
able background emission. We have applied a new method 
which, using interactively the COBRA code, is able to elim- 
inate the emission of adyacent regions and, specially, the 
low-frequency contamination. This contamination is quite 
relevant in this object due to the contribution of the un- 
derlying stellar population emission to the surface bright- 
ness of the galaxy. 

After the underlying emission was subtracted, we esti- 
mated the number of pixels above different thresholds (see 
Appendix for a complete description). From the relation 
between the number of pixels and the threshold used we 
determine the physical size of the region at the e-folding 
of the equivalent two-dimensional gaussian. The contours 
derived in this way will be not circular (see Fig. H). This 
method avoids sistematic effects introduced in the con- 
tours and sizes determination due to changes in the un- 
derlying emission from one region to another. 

The equivalent radii derived are given in Table ^ (col- 
umn 5). Finally, this procedure allows to take into account 
the effect of the Point Spread Function (PSF hereafter) 
over the sizes derived. Thus, applying = cr^oasurod ~ 
^seeing' obtaiu the radii corrected from the PSF ef- 
fect (see column 6). The e-folding radius of the PSF in 
the Ha and i?-band images was tTsooing=0'.'72. The radii 
at the e^, e'^-folding was, respectively, V2 and VS times 
the e-folding radii derived. 

4-3. Object classification 

We have spectroscopic data for only 22 of the 85 regions 
detected in the neighborhood of Mrk 86 (see Table ||) , 20 
of them in the low resolution spectra and two, #18 and 
#40, in the high resolution ones. All these regions seem to 
belong to the galaxy, with heliocentric emission-lines ve- 
locities within the range 400-540 kms~^ (see GZG). How- 
ever, the nature of the remaining 65 objects is much more 
uncertain. 

Fortunately, there are many other criteria that provide 
important clues about the nature of these objects, 1) if the 
host galaxy and some of these objects show line-emission 
(e.g., Ha emission) within the wavelength range covered 
by a narrow- band filter (FWHM~50-100A) they wiU prob- 
ably have similar recession velocities within a range of 
Au'^1000-2000kms~"'^; 2) if one of these objects is placed 
in a high surface brightness region of the galaxy, and it has 
an extended morphology, it will probably belong to the 
host galaxy. On the other hand, regions that do not show 
line emission and 3) are placed in galaxy outer regions, or 



on 




SiivX1403 
SiivX1394 CIVX1548 



1400 1600 1800 

Wavelength (A) 

Fig. 4. lUE spectrum in the range 1150-1975A obtained 
with the SWP camera in the low dispersion mode. The 
spectrum has been smoothed using a 2 pixels boxcar fil- 
ter. The SiivAA1394,1403A and CivA1548A spectral lines 
are shown with a clear P Cygni profile in the case of the 
C IVA1548A line. 



4) show point-like morphology, should not be classified as 
belonging to the host galaxy. 

In Table | we mark those regions identified spectro- 
scopically as belonging to Mrk 86 with an S letter in 
column 8. Those regions showing photometric Ha line- 
emission are marked with an E letter, and with an N those 
regions with extended morphology placed at short galacto- 
centric distances. Extended objects placed at large galac- 
tocentric distances have been marked with a B (proba- 
bly background galaxies) . Finally, three point- like sources 
have been classified as F type (probably foreground stars; 
#4, #61 & #67). 

It is worthwhile to check whether these point-like ob- 
jects without em issi on lines belong to the galaxy. Their 
colors (see Sect. 3^) indicate that objects #4, #61 and 
#67 are likely old type stars. Column 13 of Table ^ shows 
the R-band magnitudes of the regions after being sub- 
tracted from the galaxy background. Using a distance 
modulus of 29.2 we derive absolute magnitudes for these 
stars between 5 and 8 magnitudes too bright to belong to 
the galaxy. We conclude that objects #4, #61 and #67 
are effectively foreground stars. In this work we will study 
only the S, E and N type objects. 



4-.4. Optical and nIR photometry 

Using the contours obtained with COBRA, we have mea- 
sured aperture BVRJHK magnitudes and B — V, V — R, 
R — J, J — H, J ~ K colors for all the regions in Table |^ 
(columns 2-7 for magnitudes, and 8-12 for colors). These 
apertures include both knot and underlying stellar popula- 
tion emission. Before measuring these colors we degraded 



# 
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10 
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12 
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14 
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21 
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25 
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28 
29 
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31 
32 
33 
34 
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36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
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52 
53 
54 
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65 
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Table 4: Photometric magnitudes, colors and Ha fluxes. 
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Table 4: (cont.) Photometric magnitudes, colors and Ha fluxes. 
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49±0 


04 


20 


98±0 


09 


19 


89±0 


23 


19 


05±0 


17 


19 


40±0 


47 





50±0 


09 





51±0 


10 


1 


09±0 


24 





84±0 


29 


-0.35±0.50 


22 


511±0 


060 










83 


21 


46±0 


11 








20 


19±0 


11 


18 


78±0 


25 


17 


S7±0 


18 


17 


90±0 


37 














1 


41±0 


28 





91±0 


31 


-0.03±0.41 


21 


157±0 


088 










84 


22 


31±0 


11 








21 


46±0 


08 


19 


85±0 


32 


19 


00±0 


23 


19 


45±0 


71 














1 


61±0 


33 





84±0 


39 


-0.45±0.75 











305E 


15±0.277E 


17 


85 


22 


10±0 


14 








21 


14±0 


09 


19 


92±0 


60 


19 


71±0 


73 




















1 


22±0 


61 





22±0 


95 













448E 


15il:0.454E 


17 



o 
o 



^ Magnitudes and colors are not corrected for 

(1) Knot number (reverse DEC sorted) 

(2) Johnson-i? magnitude 

(3) Johnson-V magnitude 

(4) Johnson- R magnitude 

(5) J magnitude 

(6) H magnitude 

(7) 7t-prime magnitude 



extinction 



(S) B — V 

(9) V - R 

(10) R - J 

(11) J - H 

(12) H - K' 

(13) Johnson-R magnitude (underlying emission subtracted and not extinction corrected) 

(14) Ha flux (underlying emission subtracted and not extinction corrected) in ergs^^ cm 



a 
■ft 
ft 



PL, 



< 



Table 5: Emission-line fluxes. 



Knot 


Slit 


[Oil] 
A3727 


[Nelll] 
A3869 




H7 






[OIII] 
A4363 




H/3 






[OIII] 
A4959 




[OIII] 
A5007 




[NI] 
A5199 




Hel 
A5876 




[OI] 
A6300 




Ha 




[Nil] 
A6583 




Hel 
A6678 




[SII] 
A6716 




[SII] 
A6731 


E(B V) 


14 


6b 


1.25E 


15 










































































16 


3 














2 


47E 


16 


5 


27E 
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5 


07E 
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19 
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8 


21E 


16 


5 


38E 


16 


1 


35E 
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20 


4b 


3.17E 
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6 


26E 


16 








1 


OlE 


15 


2 
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20* 


4r 


1.22E 


14 
















6 


08E 


15 


3 


lOE 


15 


1 


03E 


14 






1 


31E-15 






1 


74E 


14 


3.63E 


15 






3 


30E 


15 


2 


97E 


15 


0.342 


26* 


Ir 














2 


96E 


15 


6 


36E 


14 








1 


12E 


13 






8 


78E-15 






1 


82E 


13 


2.80E 


14 


2 


94E-15 


2 


14E 


14 


1 


60E 


14 


0.277 


26* 


2r 


5.40E 


14 


3.86E-15 


1 


51E 


14 


1 


57E 


15 


3 


65E 


14 


2 


25E 


14 


6 


74E 


14 


5 


90E-16 


5 


32E-15 


9 


85E-16 


1 


04E 


13 


1.65E 


14 


1 


llE-15 


1 


16E 


14 


8 


31E 


15 


0.137 


26 


lb 


l.OOE 


13 


6.77E-16 








1 


02E 


15 


3 


80E 


14 


2 


38E 


14 




05E 


14 










































26* 


2b 


5.91E 


14 


1.99E-14 


6 


54E 


14 


2 


78E 


15 


1 


40E 


13 


7 


55E 


14 


2 


28E 


13 


1 


91E-15 




































0.534 


26* 


6b 


4.62E 


14 




3 


96E 


14 


7 


18E 


15 


8 


50E 


14 


4 


20E 


14 


1 


28E 


13 








































0.557 


27* 


6r 




















1 


94E 


14 


1 


lOE 


14 


3 


35E 


14 






1 


36E-15 


6 


41E-16 


5 


54E 


14 


9.66E 


15 






9 


02E 


15 


6 


94E 


15 


0.091 


29 


4b 


4.88E 


15 
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55E 


16 








1 


13E 


15 










































29 


3 














1 


80E 


16 


1 


70E 


15 


1 


19E 


15 


3 


55E 


15 










































32 


2b 


2.24E 


14 


6.18E-16 








2 


54E 


16 


8 


40E 


15 


4 


40E 


15 


1 


33E 


14 


2 


75E-16 






































37* 


5r 




















1 


07E 


14 


5 


57E 


15 


1 


32E 


14 










6 


35E-16 


3 


05E 


14 


6.95E 


15 






5 


89E 


15 


4 


12E 


15 


0.131 


37* 


5b 


2.53E 


14 




6 


62E 


15 








1 


42E 


14 


3 


80E 


15 


1 


12E 


14 








































0.211 


37* 


6b 


2.81E 


14 




3 


31E 


14 


1 


54E 


14 


7 


HE 


14 


2 


05E 


14 


6 


31E 
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05E 
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2r 
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14 
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6 


89E 


16 


9 


87E 


15 


5 


26E 


15 


1 


50E 


14 






1 


30E-15 


3 


78E-16 


2 


82E 


14 


5.15E 


15 






5 


18E 


15 


3 


77E 


15 


0.141 


45* 


Ir 


3.46E 


15 
















7 


43E 


15 








1 


69E 


14 














2 


12E 


14 


8.24E 


15 






7 


34E 


15 


4 


94E 


15 
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45 
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7 


45E 
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7 


67E 
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45 


lb 


5.31E 
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1 


OOE 
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45 


4b 
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7 


79E 
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1 


64E 
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80E 


15 


3 


35E 


15 


1 


09E 


14 














2 


23E 


14 
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15 






4 


41E 


15 


2 


96E 


15 


0.069 


52* 


2b 


1.92E 


14 




2 


05E 


14 


2 


OlE 


15 


4 


40E 


14 


1 


90E 


14 


5 


32E 


14 


8 
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54* 


4r 














5 


25E 


16 


1 


97E 


15 








2 


72E 


15 










4 
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5 


65E 


15 
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15 






2 


68E 


15 


2 


16E 


15 
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54 
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2 


18E 


15 
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9 


73E 


16 
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OOE 
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4b 
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2 
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1 
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2 
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62 
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16 
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16 
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Age (Myr) Age (Myr) 



Fig. 5. Model predictions. The behavior of the evolutionary synthesis models is shown. We give the change in the 
B — V and V — J colors, Ha luminosity and Ha equivalent width for two different metallicities and burst strengths. 
Solid-lines represent models with two fifths solar metallicity, and dashed-lines represent solar metallicity models. 
Thick-lines are for 6=0.01 burst strength models and thin-lines are for pure burst models. These models have been 
computed assuming a fraction of 15 per cent of escaping Lyman photons. 



the BVRHK images to the worst seeing J-band image 
(FWHM~1.8"). In addition, we measured integrated mag- 
nitudes in i?-band for these regions, subtracted from the 
underlying emission as determined by COBRA (Table ^; 
column 13). Background subtracted Ha fluxes were also 
measured and they are given in Table ^ (column 14) . Due 
to the very bad seeing of the [Oiii]A5007A image we do 
not include photometric [Oiii]A5007A fluxes in our anal- 
ysis. 

The magnitudes and colors shown in Table ^ are mea- 
sured quantities and they have not been corrected for ex- 
tinction. 



4.5. Optical spectroscopy 

The relatively small slit width employed in these obser- 
vations, prevents us from obtaining knot total emission 
line fluxes. We will derive the emission line ratios in order 
to characterize the ionized gas properties in the galaxy 
star-forming regions. In Table | we give the line intensi- 
ties measured in a region of 4'.'30x2'.'65 centered in the 
maximum of the emission knot section covered by the 
slit. These line intensities have been corrected for inter- 
nal extinction when gas color excesses, E{B — T^)gas, were 



given. The E{B — V)gaa values have been deduced from 
the Ha-H/3 and H/3-H7 Balmer line ratios measured, and 
assuming the line ratios predicted for the recombination 
case-B by Osterbrock ( 1989 ). H/3 line intensities have been 
measured deblending the absorption and emission compo- 
nents using the IRAF-STSDAS ngaussfit task. The line 
intensities were measured separately in both 6 (blue) and 
r (red) spectra. When the line-ratios measured in both 
arms, typically [O iii]A5007/H/3, matched for a given re- 
gion, we merged both data sets. 

4-.6. Ultraviolet spectroscopy 

The aperture used in the UV spectra taken from the lUE 
Final Archive (20"xl0") was centered in the brightest vi- 
sual knot (#26) with position angle PA=117°. The UV 
spectra of this knot shows the strong absorption lines of 
Si ivAA1394,1403A and C ivA1548A (see Fig. |). Although 
this spectrum has been smoothed, the signal-to-noise ra- 
tio prevents us from carrying out a quantitative analy- 
sis. However, this spectrum contains valuable information, 
showing a clear P Cygni profile on the CivA1548A line, 
typical of the fast and dense winds of 0-type supergiants. 
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Fig. 6. Color profiles in the outer regions of Mrk 86. 

5. Evolutionary synthesis models 

In order to derive the physical properties of the different 
stellar populations in Mrk 86 we have built a complete set 
of evolutionary synthesis models based in those developed 
by Bruzual & Chariot (priv. comm.). 

First, we have assumed that the optical- near-infrared 
spectral energy distribution (SED) in any region of Mrk 86 
can be described by a young star-forming burst super- 
imposed on an underlying stellar population. The burst 
strength parameter, &, will describe the mass ratio be- 
tween both stellar populations. Then, using this parame- 
ter and the evolution with time of the stellar continuum 
prediced by the Bruzual & Chariot (priv. comm.) mod- 
els, we derive the colors and number of Lyman photons 
emitted for different composite stellar populations, with b 
ranging between lO"** and 1. We have studied models with 



and 

Scalo ( |1986D IMF with Miow=O.lM0 and M^p=125Mq. 



The underlying population has been parametrized using 
the optical-near-infrared colours measured in the outer 
regions of Mrk 86. We use, B - V - i?=0.52, 

R - J=1.27, J - H=Q.m and R - K=2.'ib, as underly- 
ing stellar population colors, assuming that no significant 
color gradients are present (see Fig. In addition, the 
mass-to-light ratio adopted for this stellar population in 
the if-band was 0.87Mq/L^ q (see Paper II). 

Now, following the procedure described by Gil de Paz 
et al. (2000b) and Alonso-Herrero et al. ( 1996D we included 
the contribution to the total flux and colors arising from 
the nebular continuum and the most intense emission-lines 
(i.e., [Oii]AA3726,3729A, H/3, [Om]A4959A, [Oiii]A5007A 
and Ha, etc). 

We have assumed, in order to compute the nebular 
continuum emission, an electron density, rio, of 10^ cm^'^ 
and a temperature, Te, of 10^ K. In addition, from the 
analysis of our spectroscopic data (see Paper II), we 
adopted a A^[Heii]/A^[Hil] abundance ratio of 0.12. We 
have also assumed that the He ill abundance is so low that 
the emission from recombination to He ii is negligible. 

Finally, we have included the contribution of the emis- 
sion lines to the total flux. The contribution of the H/3, 
Ha, Pa/3, Brio-Brig and Br7 hydrogen emission lines to 
the BVRJHK bands was obtained assuming the case-B 
of recombination (Osterbrock 198£) and using the rela- 
tion given by Brocklehurst ( |1971 ) . The contribution of the 
most intense forbidden lines have been estimated using av- 
erage [Oii]AA3726,3729/[Oiii]A5007 and [O iii]A5007/H/3 
line ratios, as provided by our spectroscopic data. Fortu- 
nately, the contribution of all the forbidden lines to the 
B and V bands is very small. Using the higher and lower 
line-ratios measured in the galaxy, this contribution would 
range between 1 and 8 per cent for the 5-band and 2 and 
8 per cent for the F-band, for a Ha equivalent width (EW 
hereafter) of lOOA. 

The output of the models will be the optical-near- 
infrared colors B~V,V~R,V-J,J~HandJ~Kof 
the composite stellar population, its Ha luminosity and 
equivalent width, and mass-to-light ratio, parametrized 
as a function of the burst age, burst strength and stel- 
lar metallicity {t,b,Z). 

The Cousins-i? magnitudes originally given by the 
Bruzual & Chariot (priv. comm.) models have been con- 
verted to Johnson- i? magnitudes using the relation given 



by Fernie (1983). However, if we compare the correction 
predicted by Fernie ( 1983| ) in the case of very red stars 
{Rc ~ i?j~0.25™) with that measured by Fukugita et al. 
( 1995[ ) for early-type galaxies, typically of 0.1™ -with no 
correction for extinction applied-, we find differences of 
about 0.15™. Since the change in Rc — Rj due to the cor- 
rection for extinction can not be higher than 0.02™, this 
difference has to be attributed to a difference in the cor- 
rection between evolved stellar populations and individual 
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very-red stars. Thus, in the case of the underlying popu- 
lation analysis, we have applied the mean correction given 
by Fukugita et al. ( |1995| ). 



6. Results 

6.1. Underlying stellar population 

The surface brightness profiles given in Fig. || show a clear 
exponential dominance at regions outer than '-^1.25kpc. 
We have obtained the parameters of this exponential com- 
ponent on the deepest i?-band profile yielding a scale, a, of 
0.88±0.02kpc and an extrapolated central surface bright- 
ness, /iE.Oj of 21.50±0.06magarcsec~^. Papaderos et al. 
( 1996a ) obtained quite different values, a=0.74±0.06kpc 
and /iE,o=19.42±0.03magarcsec~^. This difference is not 
surprising considering that the i?-band profile given Pa- 
paderos et al. (1996a) only reaches galactocentric dis- 
tances of 2 kpc. 

The color profiles in the galaxy outer regions, where 
the underlying stellar population dominates the total light 
profile, are shown in Fig. |[ No significant color gradients 
are observed at distances larger than 1.2-1.3kpc, whereas 
at shorter distances a progressive and swallow blueing is 
derived. This blueing is probably related with an incre- 
ment in the light contamination from high surface bright- 
ness star-forming regions associated to the plateau com- 
ponent. In Fig. we show the Ha profile compared with 
the surface brightness profile in the S-band of the plateau 
component (see Papaderos et al. 1996a). From this figure 
it is quite clear that the progresive blueing observed at dis- 
tances shorter than 40" is due to currently star-forming 
regions located in the plateau component. 



6.2. Star-forming regions 



In Sect. 4.3 we have classified the regions observed in the 
Mrk 86 neighborhood as S, E, N, F and B objects (S, 
spectroscopically confirmed; E, emission-line regions; N, 
extended and diffuse; F, foreground stars; B background 
galaxies) . 

The optical-nIR colors measured for the F objects 
seem to indicate that the #4, #61 and #67 regions are, 
respectively, M0-M3, K2-K5 and G7-K0 spectral type fore- 
ground stars. Besides, between the B type objects (back- 
ground galaxies), we find two very red objects, #2 and #3 
regions , with V—K colors of about 5 and 4™, respectively. 

Hereafter, we present the results for the study of the 
S, E and N-type regions. First, we have analyzed their 
i?-band luminosities in Fig. |8^. We observe that this dis- 
tribution has a clear maximum at M/j=— 9.5™. This dis- 
tribution can not be well fitted using any standard power- 
law luminosity function (see, e.g. Elson & Fall |1985| , for 
the LMC star-clusters LF). 

In Fig. [sja we show the Ha luminosity distribution of 
the galaxy emission-line regions. The Ha luminosities used 




20 40 60 80 
Radial distance (arcsec) 



100 



Fig. 7. Ha profile obtained using the isophotes measured 
in the i?-band image {solid-line). The surface brightness 
profile of the plateau component as parametrized by Pa- 
paderos et al. ( 1996a ) in the i?-band is also shown [dot- 
dashed-line) . 



were corrected for internal extinction. In those cases were 
spectroscopic data were not available we used a mean color 
excess of E{B — y)gas=0.34™. The distribution obtained 
is nearly flat in the luminosity range lO'^^'^-lO^^'^ergs"^, 
with most of the H ii regions showing Ha luminosities 
in the range lO'^^-lO'^^ ergs"^. This luminosity range cor- 
responds to star formation rates between 6xl0~^ and 
O.O6M0yr~i (Gil de Paz et al. pOOObj ). The spatial res- 
olution of the Ha image is about 40 pc. Therefore, the 
scarcity of H ii regions fainter than 10'^^ ergs""'^ can not be 
explained as an e ffect of the spatial resolution (see Ken- 
nicutt et al. 



It should be noticed that some of these Ha emitting 
regions showing very faint or null continuum emission (#7, 
#9, #12, #50), seem to be pure gas regions photoionized 
by distant stellar clusters. These regions are in many cases 
associated with the rim of expanding bubbles (see Martin 



19981 , GZG). 

Now, in Fig. ^ we show the radius distribution mea- 
sured in the i?-band image (given by the e-folding radius; 
see Appendix). These radii have been corrected for the 
PSF contribution (see Sect. This distribution shows 
a maximum at about 1 arc second, that corresponds to a 
FWHM of about 55 pc. 

In Fig. ^d, the i?-band and Ha luminosities are com- 
pared. The lines drawn suggest that these regions have 
equivalent widths of Ha ranging between 100-500A. bmce 
both. Ha and i?-band fluxes have been measured sub- 
tracted from the contribution of the underlying popula- 
tion, these EW(Ha) are pure H II region equivalent widths. 
In this figure the symbol size is proportional to the extinc- 
tion corrected B — V color, using larger symbols for bluer 
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Table 6. Data for the Mrk 86 starburst component. 
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Fig. 8. a) i?-band luminosity distribution for the S, E 
and N-type regions, b) Ha luminosity distribution for the 
galaxy emission- line regions (S and E types), c) Radius 
distribution for the S, E and N-type regions, d) Ha and 
continuum i?-band luminosities comparison. Lines of con- 
stant EW(Ha) are drawn, e) i?-band apparent magnitude 
vs. radius. 



Coordinates 


RA(J2000) 


08" 13" 14.69" 


DEC(J2000) 


-(-45°59'21.9" 


Photometry 


B-V 


0.39±0.06 


V -R 


0.54±0.25 


V - J 


1.59±0.02 


J-H 


0.77±0.08 


J-K' 


0.85±0.03 


Mr 


-14.0±0.2 


Radius 


4.3" 




Spectroscopic indexes 


SUt 


#4b #lb 


D4000 


1.38 1.38 


Mg2 


0.06 0.04 




6.0A 6.2A 


Fe5270 


I.20A 


Fe5406 


O.74A 


Region 


21'.'45x2'.'65 14'.'3x2'.'65 



^ Equivalent width in absorption. 



regions. Finally, we have compared (see Fig. ||e) the knot 
apparent magnitude and the physical radius, both mea- 
sured on the i?-band image using the COBRA program. If 
all these star-forming regions were optically thin at these 
wavelengths and they had similar star densities we would 
expect the flux to be proportional to the cube of the ra- 



dius. In Fig. 
also shown. 



the best fit to the function F» (x B? is 



6.3. Central Starburst 

Finall y, we h ave studied the region classified by Papaderos 
et al. ( 1996a ) as the starburst component (see Fig. ^). This 
component appears as a bright, extended and not very well 
defined region in the broad-band BVRJHK images (see 
Fig. |l|). 

The integrated blue color measured by Papaderos et al. 
( |1996aD in the galaxy central parts, B-R ~0.9 (0.93±0.26 
in our work) , similar to that observed in Scd and Im galaxy 
types (Fukugita et al. 1995), indicates the existence of a 
young stellar population superimposed on the evolved un- 
derlying component. However, the absence of significant 
gas emission (see the Ha image in Fig. |l]) suggests an 
intermediate aged dominating population. In Table ^ we 
give the colors and spectr oscop ic indexes measured for this 
region (see Trager et al. 1998 ) at two different slit posi- 
tions^. The integrated optical-near-infrared colors shown 
in Table |^ have been measured on the i?-band image using 
an aperture with radius 4.3" (1 e-folding), as given by the 
COBRA program. This aperture is marked in Fig. with a 



^ The spectroscopic indexes hav e bee n measured using the 
REDUCEME package (Cardiel et al. 1998). 
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thick-lined contour. These colors have not been corrected 
for extinction. The Mr absolute magnitude has been mea- 
sured subtracted from underlying emission. The spectro- 
scopic indexes have been measured in the integrated spec- 
tra of the starburst section covered by the slits #4b and 
#lb. (see Table |; see also Fig. 1 of GZG). 

In Paper II we will derive the physical properties of 
this region by comparing the optical-nIR colors obtained 
with the predictions of our evolutionary synthesis models. 
In addition, the spectroscopic indexes measured will be 
compared with those prediced by the Bruzual & Chariot 
(priv. comm.) evolutionary synthesis models. 

7. Conclusions 

In this paper, the first of two, we have presented the ob- 
servations and data analysis for the optical-near-infrared 
study of the Blue Compact Dwarf galaxy Mrk 86. We have 
taken BVRJHK, Ha and [Oiii]A5007A images and long- 
slit optical spectroscopy for the object. Thus, using all 
these data, 

— we have obtained the underlying population surface 
brightness and color profiles. The deepest i?-band 
image yields an exponential profile for this compo- 
nent with extrapolated central surface brigthness of 
ME,o=21.5magarcsec~^ and scale, a=0.88kpc. No sig- 
nificant color gradients are observed in this component 
at distances larger 1.2-1.3kpc. 

— We have cataloged and classified all the regions ob- 
served in the neighborhood of the galaxy. If these re- 
gions are associated with the galaxy have been classi- 
fied as S (spectroscopically confirmed), E (confirmed 
by their photometric Ha or [Oiii]A5007A emission) or 
N (diffuse and placed in the galaxy central region) type 
regions. S and E type objects are accepted to be galac- 
tic Hii regions. N objects are, probably, evolved stel- 
lar clusters, with negligible Lyman continuum emis- 
sion. On the other hand, background galaxies and fore- 
ground stars are classified as B and F types, respec- 
tively. 

— We have also described the set of evolutionary syn- 
thesis models used in Paper II. These models are 
based in those developed by Bruzual & Chariot (priv. 
comm.). They have been obtained for metallicities be- 
tween 1/50 Zq and 2Zq, and burst strength ranging 
between 10~^ and 1. We have included nebular contin- 
uum and recombination and forbidden lines emission. 
Although the contribution of the forbidden-line fluxes 
is quite uncertain in those regions where spectroscopic 
data are not available, the contribution of these lines 
in the RJHK bands is negligible and lower than 8 per 
cent for the i?T^-bands (assuming an EW(Ha)~100A). 

— The optical-near-infrared colors, fluxes, emission line 
intensities of the S, E and N regions have been mea- 
sured. Optical and near-infrared colors have been mea- 




10 lb 
pixels 

Fig. 9. Example of an observed x-axis light profile {thick 
solid-line) and its corresponding best fit (thin solid-line). 
The dotted-lines correspond to the two gaussian compo- 
nents, being the dot-dashed-line the sum of these two com- 
ponents. The underlying emission (background) is repre- 
sented by the dashed-line. The approximate positions for 
the 8 points used for defining the initial set of parameters 
are also shown. The spatial scale of the image was 0'.'333 
pixel" ^. 



sured. The physical properties of the galaxy star- 
forming regions will be derived in Paper II, compar- 
ing these data with the evolutionary synthesis models 
developed. This comparison will be performed using 
a method based in the combination of Monte Carlo 
simulations, a maximum likelihood estimator, a single 
linkage clustering anaylisis method and the Principal 
Component Analysis (PCA). 
— Finally, we provide optical-near-infrared colors and 
spectroscopic indexes for the central starburst com- 
ponent. 

8. Appendix: COBRA program 

Due to the large number of star-forming regions in Mrk 86 
and the intense and variable underlying emission we de- 
veloped a program that allowed us to substract this 
underlying emission and the occasional contamination 
from neighbor star-forming regions. Then, on the un- 
derlying emission subtracted image, we determined the 
most reliable apertures, sizes, and integrated fluxes for all 
these regions. This program was written in FORTRAN?? 
programming language and it was called COBRA (see 
tittp : //www . ucm . es/ info/ Astrof /cobra/) . 

COBRA uses different graphic output devices, includ- 
ing xwiNDOWS and Pos tscrip t, and makes use of the 
AMOEBA (Press et al. |l986D , PGPLOT, F ITSIO and 
BUTTON (Cardiel et al. |1998| see also Cardiel |l999| ) sub- 
routines. 

In this Appendix we will describe the procedure fol- 
lowed to estimate the underlying emission and the criteria 
used to determine the sizes and fluxes of the star-forming 
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regions and stellar clusters in Mrk 86 using the program 
COBRA. 



8.1. Light profiles fitting 

First, (1) a image section including the region of interest is 
selected interactively (see the original images in Fig. |lO| ). 
The size of this section necessarily depends on the spatial 
scale of the change in the underlying emission and the 
proximity of neighbor regions. This section must be chosen 
in such a way that the background emission (underlying 
and neighbor regions emission) can be reproduced by a line 
in the x and y light profiles. Thus, when several relative 
maxima are placed too close to match this criterion they 
were studied as a single region. 

After cropping the image section selected, (2) the po- 
sition of the maximum of the region of interest is marked 
interactively. Then, the light profile along the x axis at 
this position is fitted. 

In order to perform this fit three different components 
are used, a line that reproduced the underlying and neigh- 
bor regions contamination and two gaussians that allowed 
to reproduce the light profile of the region of interest. In 
principle, any smooth, monotonic light profile can be ap- 
proximated by a serie of Gaussians (somewhat analogous 
to the Fourier series), not only the PSF of an image, but 
also the light profiles of remote stellar clusters or galac- 

199(]| ). Two gaussians are 



tic cores (see Bendinelli et al. 

employed in order to reproduce the light profile of these 
regions because the use of a larger number of components 
led in some cases to solutions with no physical meaning. 

Then, the best fit for the initial light profile using 
the minimization subroutine AMOEBA is obtained. The 
AMOEBA subroutine needs an initial set of solutions 
which can be defined by the 8 parameters of the func- 
tion to minimize (3 for each gaussian and 2 for the line). 
These parameters were introduced interactively marking 
the approximate position of the center and FWHM for 
each gaussian and two extreme values in the profile for 
the line. In Fig. || we show an example of the positions of 
the points used to define the set of initial parameters for 
the fit. Points 1,2,3 and 4,5,6 led to the initial parameters 
for the two gaussians and points 7,8 for the line. Using 
this set of parameters, (3) the best-fitting light profile is 
obtained (see Fig. ||). 

Then, from the output of this fit and introducing the 
new positions of the two gaussian maxima, (4) all the x 
axis profiles are fitted starting from those profiles adjacent 
to the initial one. The need of using previous results for 
the subsequent fits constitutes the main reason for starting 
this fitting procedure at the maximum of the emission 
region. 

Next, (5) the same fitting procedure is applied to the y 
axis profiles. Finally, after the profiles have been fitted in 
both axis, the underlying emission images reconstructed 
are averaged and subtracted from the original input image 
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Fig. 11. Change in the area {A^) of the region #18 with 
the threshold, The solid-line represents the result ex- 
pected for a radius of 0'.'9 at 1 e-folding. The dashed-lines 
correspond to 0'.'8 and I'.'O e-folding radii. 



(see several examples in Fig. |i^). This background- free 
image is then used to determine the centers, apertures, 
sizes and integrated fluxes of the different regions. 

8.2. Position, sizes and fluxes 

The knot center is taken as the maximum of the sum of 
the two gaussian components. Then, computing the knot 
surface, A^, above different thresholds, /x, we can estimate 
the equivalent knot e-folding radius, i?e-foiding: using. 



-^c — folding 



7rln(i^) 



(1) 



where /max is the knot maximum intensity. In Fig. ^ we 
show the change in the area A^ for different thresholds as 
measured in the background subtracted image of the star- 
forming region #18 (with 1 e-folding radius of 0'.'93). The 
relatively small change deduced for the e-folding radius 
with the threshold guarantees the goodness of this size 
determination procedure. 

The e^ and e^-folding radii are then computed as a/2 
and -v/S times the e-folding radius. Then, the knot contours 
are derived at these three radii, i.e. at /niax//x=e,e^,e'^. 
The knot total fiux is computed as the total flux in the 
background subtracted image. In order to ensure that 
most of the knot fiux was included we derive the knot 
growing-curve, finding differences not larger than 20 per 
cent between total and growing-curve extrapolated flux. 
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Fig. 10. Original, total (reconstructed), background, original minus background model and residuals images in both 
axis and the average as obtained for the region #6 in the i?-band. The cuts in the original, total and background 
images are the same and equal to the maximum and minimum value in the original image. The cuts in the residuals 
images correspond to plus/minus one tenth the difference between the maximum and minimum in the background 
averaged image. The residuals images were obtained as the original image minus the total (reconstructed) image. 
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Fig. 10. (cont.) Original, total (reconstructed), background, original minus background model and residuals images in 
both axis and the average as obtained for the region #11 in the -R-band. The cuts in the original, total and background 
images are the same and equal to the maximum and minimum value in the original image. The cuts in the residuals 
images correspond to plus/minus one tenth the difference between the maximum and minimum in the background 
averaged image. The residuals images were obtained as the original image minus the total (reconstructed) image. 
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Fig. 10. (cont.) Original, total (reconstructed), background, original minus background model and residuals images in 
both axis and the average as obtained for the region #18 in the -R-band. The cuts in the original, total and background 
images are the same and equal to the maximum and minimum value in the original image. The cuts in the residuals 
images correspond to plus/minus one tenth the difference between the maximum and minimum in the background 
averaged image. The residuals images were obtained as the original image minus the total (reconstructed) image. 
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Fig. 10. (cont.) Original, total (reconstructed), background, original minus background model and residuals images in 
both axis and the average as obtained for the region #70 in the -R-band. The cuts in the original, total and background 
images are the same and equal to the maximum and minimum value in the original image. The cuts in the residuals 
images correspond to plus/minus one tenth the difference between the maximum and minimum in the background 
averaged image. The residuals images were obtained as the original image minus the total (reconstructed) image. 



A. Gil de Paz et al.: Mapping the star formation history of Mrk 86* 



21 



Ultraviolet Explorer archive at the ESA VILSPA observa- 
tory. Based also on observations collected at the German- 
Spanish Astronomical Center, Calar Alto, Spain, operated 
by the Max-Planck- Institut fiir Astronomic (MPIA) , Hei- 
delberg, jointly with the Spanish 'Comision Nacional dc 
Astronomia'. 

We arc grateful to Carmc Jordi and D. Galadi for ob- 
taining the F-band image. We would like to thank C. 
Sanchez Contreras and L.F. Miranda for obtaining the 
high resolution spectra, N. Cardiel for providing the re- 
DUCEME package and C.E. Garcia Dabo for the PlateAs- 
TROM routine. We also thank A. Alonso-Herrero for her 
help in the acquisition and reduction of the nIR images. 
We also acknowledge the anonymous referee for several 
helpful comments. Finally, we thank to C.E. Garcia Dabo 
and J. Cenarro for stimulating conversations, and M. Sha- 
rina for providing a reprint copy of her article. This re- 
search has been supported in part by the grants PB93-456 
and PB96-0610 from the Spanish 'Programa Sectorial de 
Promocion del Conocimiento'. A. Gil de Paz acknowledges 
the receipt of a 'Formacion del Profesorado Universitario' 
fellowship from the Spanish 'Ministerio de Educacion y 
Cultura'. 

References 

Alonso-Herrero A., Aragon-Salamanca A., Zamorano J., Rego 

M., 1996, MNRAS 278, 417 
Aragon-Salamanca A., Ellis R.S., Couch W.J., Carter D., 1993, 

MNRAS 262, 764 
Arp H., 1966, Atlas of Peculiar Galaxies. California Institute 

of Technology, Pasadena 
Bendinelli O., Paxmeggiani G., Zavatti F., Djorgovski S., 1990, 

AJ 99, 774 

Bonatto C., Bica E., Pastoriza M.G., AUoin D., 1999, A&A 

343, 100 

BottincUi L., Gouguenheim L., Fouque P., Paturel G., 1990, 

A&AS 82, 391 
Brocklehurst M., 1971, MNRAS 153, 471 
Burkcrt A., 1989, Ph.D. Thesis, Munich 
Campbell A.W., Terlevich R., 1984, MNRAS 211, 15 
Cardiel N., 1999, Ph.D. Thesis, Universidad Complutense de 

Madrid 

Cardiel N., Gorgas J., Cenarro J., Gonzalez J. J., 1998, A&AS 
127, 597 

Cassatella A., Barbero J., Benvenuti P., 1985, A&A 144, 335 
Doublier V., 1998, Ph.D. Thesis, Marsella 

Drinkwater M., Hardy E., 1991, AJ 101, 94 
Dultzin-Hacyan D., Masegosa J., Moles M., 1990, A&A 238, 
28 

Elias J.H., Frogel J.A., Matthews K., Neugebauer G., 1982, AJ 
87, 1029 

Elson R.A.W., Fall S.M., 1985, PASP 97, 692 
FaneUi M.N., O'ConcU R.W., Thuan T.X., 1988, ApJ 334, 665 
Fernie J.D., 1983, PASP 95, 782 
Fuentes-Masip O., 1997, Ph.D. Thesis, La Laguna 
Fukuguita M., Shimasaku K., Ichikawa T., 1995, PASP 107, 
945 

Garci-Dabo C.E., Gallego J., 1999, Ap&SS 263, 1-4 



Gil de Paz A., Zamorano J., Gallego J., 1999, MNRAS 306, 

975 (GZG) 

Gil de Paz A., Zamorano J., Gallego J., 2000a, A&A (Paper II) 
Gil de Paz A., Aragon-Salamanca A., Gallego J., et al., 2000b, 

MNRAS 316, 357 
GiovancUi R., Haynes M.P., Herter T., et al., 1997, AJ 113, 53 
Hoffman G.L., Salpctcr E.E., Hclou G., 1990, in the proceed- 
ings of the Edwin Hubble Centennial Symposium, 67 
Kennicutt R.C., Edgar B.K., Hodge P.W., 1989, ApJ 337, 761 
Kent S.M., 1985, PASP 97, 165 

Klein U., Wielebinski R., Thuan T.X., 1984, A&A 141, 241 
Klein U., Weiland H., Brinks E., 1991, A&A 246, 323 
Kunth D., Sargent W.L.W., 1986, ApJ 300, 496 
Kunth D., Maurogordato S., Vigroux L., 1988, A&A 204, 10 
Landolt A.U., 1973, AJ 78, 959 

Loose H.-H., Thuan T.X., 1985, The Morphology and Structure 
of BCDGs from CCD Observations, in: Kunth, D., Thuan, 
T.X., and Van, J.T.T.(eds.) Star-Forming Dwarf Galaxies. 

Editions Fronticrcs. 
Markarian B.E., 1969, Astrofizika 5, 443 

Marlowe A.T., Heckman T.M., Wyse R.F.G., Schommer R., 

1995, ApJ 438, 563 
Martin C. L., 1998, ApJ 506, 222 
MazzarcUa J.M., Boroson T., 1993, ApJS 85, 27 
Norton S., Salzer J. J., 1997, BAAS 190 

Osterbrock D.E., 1989, Astrophysics of Gaseous Nebulae and 
Active Galactic Nuclei. University Science Books, Mill Val- 
ley, California. 

Papaderos P., Loose H.-H., Thuan T.X., Fricke K.J., 1996a, 

A&AS 120, 207 
Papaderos P., Loose H.-H., Fricke K.J., Thuan T.X., 1996b, 

A&A 314,59 

Petrosian A.R., Boulesteix J., Comte G., Kunth D., LeCoarer 

E., 1997, A&A 318, 390 
Press W. H., Flannery B.P., Teukolsky S.A., Vetterling W.T., 

1986, Numerical Recipes, Cambridge University Press 
Rood H.J., Dickel J.R., 1976, ApJ 205, 346 
Sage L.J., Salzer J.J., Loose H.-H., Henkel C, 1992, A&A 265, 

19 

Scalo J.M., 1986, Fund. Cosmic Phys. 11, 1 
Searle L., Sargent W.L.W., 1972, ApJ 173, 25 
Searle L., Sargent W.L.W., Bagnuolo W.G., 1973, ApJ 179, 
427 

Shapley H., Ames A., 1932, Ann. Harvard College Obs. 88, No. 
2 

Sharina M.E., Karachentsev I.D., Tikhonov N.A., 1999, As- 
tronomy Letters 25, 322 
Silk J., Wyse R.F.G., Shields G.A., 1987, ApJ 322, L59 
Thronson Jr. H.A., Telesco CM., 1986, ApJ 311, 98 
Thuan T.X., 1983, ApJ 268, 667 

Thuan T.X., 1991, Observations and Models of Blue Compact 
Dwarf Galaxies, in: Lcithcrcr C, Walborn R.N., Heckman 
T.M., Norman C.A. (cds.) Massive Stars in Starbursts, 
Cambridge University Press, p. 183. 

Thuan T.X., Martin G.E., 1981, ApJ 247, 823 

Trager S.C., Worthey G., Faber S.M., Burstein D., Gonzalez 
J.J., 1998, ApJS 116, 1 

Verter F., 1985, ApJS 57, 261 

Wainscoat R.J., Cowie L.L., 1992, AJ 103, 332 



This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0007267vl 



This figure "fig2.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0007267vl 



